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^ • Abstract. We show how to obtain the probability density fimction for the amplitude 

jy^ I of the curvature perturbation, TZ, produced during an epoch of slow-roll, single-field 

^ ' inflation, working directly from n-point correlation functions of TZ. These n-point 

functions are the usual output of quantum field theory calculations, and as a result 

K^ ^ we bypass approximate statistical arguments based on the central limit theorem. Our 

^ • method can be extended to deal with arbitrary forms of non-Gaussianity, appearing 



C^ 



at any order in the n-point hierarchy. We compute the probability density for the 
total smoothed perturbation within a Hubble volume, e, and for the spectrum of e. 
When only the two-point function is retained, exact Gaussian statistics are recovered. 
When the three-point function is taken into account, we compute explicitly the leading 
slow- roll correction to the Gaussian result. 
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1. Introduction 

There is now a good deal of observational evidence that the generic predictions of the 
inflationary scenario are realised in the spectrum of density perturbations in our universe 
^|21IH1I1]- For slow-roll inflation driven by a scalar 0, these predictions are: 

(i) a nearly scale-invariant spectrum of fluctuations on all scales accessible to 
cosmological observation; 

(ii) for inflation near the theoretically motivated energy scale Mqut — 10^^ GeV, these 

fluctuations should have magnitude 5p/ p ^ 10~^; and 
(iii) the fluctuation spectrum should exhibit Gaussian statistics, in the sense that the 

probability distribution of the density fluctuation should be approximately normally 

distributed. 

(For a review of the inflationary paradigm and its predictions, see eg., Refs. |Sl Ej-) 
Predictions (ji)) and (Jul) can be obtained using the standard techniques of quantum field 
theory. This calculation is now classical [3 E] and relies only on the fact that the 
vacuum fluctuation of a scalar field in de Sitter space with Hubble parameter H is 
roughly dcj) = H/27i jHl 13 CDl IH] ■ The third prediction — that the fluctuation spectrum 
is Gaussian — is less transparent. It follows from the fact that the inflaton perturbation, 
which is commonly expressed in terms of the so-called comoving curvature perturbation^ 
[T3] TZ, is treated as a free field, 

7^(t,x)= /-^7^(t,k)e^'^■^ (2) 



J (27r)3 

where there is no coupling between the TZ{t, k) for different k. With this understanding, 
Eq. (j21) means that 7?. does not interact either with itself, or any other particle 
species in the universe. The real-space field TZ{t, x) is obtained by summing an 
infinite number of independent, identically distributed, uncorrelated oscillators. Under 
these circumstances the Gaussianity of TZ{t, x) follows from the central limit theorem 
[Hj, given reasonable assumptions about the individual distributions of the TZ(t,k). 
The exact form of the distributions of the TZ{t, k) is mostly irrelevant when making 
statements about the inflationary density perturbation. 

I TZ{t, x) expresses the relative expansion of a given local neighbourhood of the universe with respect 
to the background, in a gauge where observers in free-fall with the expansion see no net momentum 
flux. This is the so-called comoving gauge. In particular, during inflation the scalar field fluctuation is 
zero on comoving slices. The metric in this gauge is thus 

ds^ ^ ~N^{t,x)dt^ + a^{t)e^'^^'^'''>-fij{dx' + N\t,x)dt){dx^ + N^{t,x)dt), (1) 

where 7^ is the metric on unperturbed spatial slices, which we take to be flat and is subject to the 
condition det7 = 1; and a{t) is the unperturbed scale factor of the universe. The functions A^ and iV' 
are the so-called lapse function and shift vector, and are determined by algebraic constraint equations 
in terms of the matter content and the metric fields {TZ, a and 7), once the gauge is fixed. Our 
sign convention for the metric is "mostly plus", (—,+,+,+). The advantage of working with the 
perturbation 7?.(i,x) is that it mixes scalar fluctuations from the metric and matter sectors in a gauge 
invariant way ^12j. 
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In conventional quantum field theory, all details of TZ and its interactions are 
encoded in the n-point correlation functions of 7^ (or their Fourier transforms), written 
(out|7^(ti, xi) ■ ■ ■7^(t„,x„)|in). Working in the Heisenberg picture, where the fields 
carry time dependence but the states {|in), |out)} do not, these functions express 
the amplitude for the early-time vacuum |in) to evolve into the late-time vacuum 
|out) in the presence of the fields 7^(tj,Xj). Given the n-point functions for all n at 
arbitrary x and t, one can determine TZ{t,x) [15 , at least in scattering theory. In the 
context of the inflationary density perturbation, these vacuum evolution amplitudes 
are not directly relevant. Instead, one is interested in the equal time expectation values 
(in|7^(t, xi) ■ ■ ■ 7l(t, x„)|in), which can be used to measure gravitational particle creation 
out of the time-independent early vacuum |in) during inflation. These expectation values 
are calculated using the so-called "closed time path formalism," which was introduced 
by Schwinger [J^ (see also [13 CHI UHl 120] ) • In this formalism there is a doubling of 
degrees of freedom, which is also manifest in finite temperature calculations [211 1^ . 
This method has recently been used |2Sl 1211 12^] to extend the computation of the 
correlation functions of TZ beyond tree-level. 

Knowledge of the expectation values of 7?. in the state |in) is sufficient to predict a 
large number of cosmological observables, including the power spectrum of the density 
perturbation generated during inflation [HI UHl , and the two- and three-point functions 
of the cosmic microwave background (CMB) [2nil2Zll2Sll2nilSIIlinillS21ISSllSllSSllSE] 
temperature anisotropy. Because they are defined as expectation values in the quantum 
vacuum these observables all have the interpretation of ensemble averages, as will be 
discussed in more detail below. On the other hand, it is sometimes required to know 
the probability that fluctuations of some given magnitude occur in the perturbation TZ 
[SZl dH ISHl • This is not a question about ensemble averages, but instead is concerned 
with the probability measure on the ensemble itself. As a result, such information cannot 
easily be obtained from inspection or simple manipulation of the n-point functions. 

For example, if we know by some a priori means that 7?- is free, then the argument 

given above based on the central limit theorem implies that at any position x, the 

probability of fluctuations in 7?. of size e must be 

1 / e^ \ 

P(7?. has fluctuations of size e) ~ — exp I —77^ , (3) 



where the variance in 7^ is 

a^ = (7^(t,x)2)= / d\n kV{k). (4) 



The quantity T^{k) is the so-called dimensionless power spectrum, which is deflned in 
terms of the two-point function of 7^, calculated from the quantum fleld theory in- 
vacuum: 

(in\TZ{tMMtM)H = {2nfS{ki + k,)'^V{k,). (5) 

This is the only relevant observable, because it is a standard property of free flelds 
that all other non- vanishing correlation functions can be expressed in terms of the two- 
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point function Q, and hence the power spectrum. In practice, in order to give a 
precise meaning to Q, it would be necessary to specify what it means for TZ to develop 
fluctuations of size e, and whether it is the fluctuations in the microphysical field TZ or 
some smoothed field TZ which are measured. These details affect the exact expression 
(^ for the variance of e. 

The average in Eq. (0), denoted by (in|---|in), is the expectation value in 
the quantum in-vacuum. To relate this abstract expectation value to real-world 
measurement probabilities, one introduces a notional ensemble of possible universes, 
of which the present universe and the density fluctuation that we observe is only one 
possible realization (see, eg., inS])- However, for ergodic processes, we may freely trade 
ensemble averages for volume averages. If we make the common supposition that the 
inflationary density perturbation is indeed ergodic, then we expect the volume average 
of the density fluctuation to behave like the ensemble average: the universe may contain 
regions where the fluctuation is atypical, but with high probability most regions contain 
fluctuations of root mean square (rms) amplitude close to a. Therefore the probability 
distribution on the ensemble, which is encoded in Eq. (0), translates to a probability 
distribution on smoothed regions of order the horizon size within our own universe. 

In order to apply the above analysis, it was necessary to know in advance that TZ 
was a free field. This knowledge allowed us to use the central limit theorem to connect 
the correlation functions of TZ with the probability distribution (jH)). The situation in 
the real universe is not so simple. In particular, the assumption that during inflation 
TZ behaves as a free field, and therefore that the oscillators 7?-(k) are uncorrelated and 
independently distributed, is only approximately correct. In fact, TZ is subject to self- 
interactions and interactions with the other constituents of the universe, which mix 
k-modes. Consequently, the oscillators 7^(k) acquire some phase correlation and are 
no longer independently distributed. In this situation the central limit theorem gives 
only approximate information concerning the probability distribution of 7?.(x), and it is 
necessary to use a different method to connect the correlation functions of TZ with its 
probability distribution. 

In this paper we give a new derivation of the probability distribution of the 
amplitude of fluctuations in TZ which directly connects P(e) and the correlation functions 
(7?.(ki) ■ ■ ■7?.(k„)), without intermediate steps which invoke the central limit theorem 
or other statistical results. When the inflaton is treated as a free field, our method 
reproduces the familiar prediction of Gaussian statistics. When the inflaton is not 
treated as a free field, the very significant advantage of our technique is that it is possible 
to directly calculate the corrections to P(e). Specifically, the interactions of TZ can be 
measured by the departure of the correlation functions from the form they would take if 
TZ were free. Therefore, the first corrections to the free-field approximation are contained 
in the three-point function, which is exactly zero when there are no interactions. 
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The three-point function for single-field, slow-roll inflation has been calculated by 
Maldacena [IHl, whose result can be expressed in the form |1T] 

(7^(kl)7^(k2)7^(k3)) = A7i\27,f5{Y^]^,)^A{hMM), (6) 

where A is one-half§ Maldacena's ^-function jlU] and V^ (to be defined later) measures 
the amplitude of the spectrum when the kj crossed the horizon. (For earlier work, see 
Refs. jini EH ESI HE] • The present situation is reviewed in Ref. jlZj.) This result has 
since been extended to cover the non-Gaussianity produced during slow-roll inflation 
with an arbitrary number of fields |iniEIlEHlEniEniEIlE2lES!, preheating PI ESI EH] , 
models where the dominant non-Gaussianity is produced by a light scalar which is a 
spectator during inflation |SZ1 EHl tHH] , and alternative models, often involving a small 
speed of sound for the inflaton perturbation |121 EHl EOl EB ESj • For single-field, slow-roll 
inflation, the self-interactions of TZ are suppressed by powers of the slow-roll parameters. 
This means that the correction to Gaussian statistics is not large. In terms of the A- 
parametrized three-point function ©, this is most commonly expressed by writing 

5 A 
/nl = -- y^ ,3 = 0(slow roll), (7) 

where /nl |2H1 E2] expresses the relative contribution of a non-Gaussian piece in 7?., 
viz, TZ = TZg — |/nl * TZ?g, and TZg is a Gaussian random field. (There are differing 
sign conventions for /nl ISHI-) In models with more degrees of freedom it is expected 
to be possible to obtain very much larger non-Gaussianities, perhaps with /nl ~ 10 
[68| EH ESI ESI EH EOl EH] , although no unambiguous concrete example yet exists except 
where the non-Gaussianity is generated during preheating jSHESj- If the inflationary 
perturbation has a speed of sound different from unity then large non-Gaussianities 
may also appear |121, although in this case it is difficult to simultaneously achieve scale 
invariance. The current observational constraint can approximately be expressed as 
|/nl| ^ 100 T]. In the absence of a detection, the forthcoming Planck Explorer mission 
may tighten this constraint to |/nl| ^ 3 pHllHlj . 

Non-Gaussian probability distributions have been studied previously by several 
authors. The closest to the method developed in this paper include that of Matarrese, 
Verde & Riotto jSZI, who worked with a path integral expression for the density 
fluctuation smoothed on a scale R (which they denoted '5r'); and that of Bernardeau 
& Uzan ISHlEni- The latter analysis has some features in common with our own, being 
based on the cumulant generating function, and moreover since the expression for the 
probability density in Refs. [SHI EH] is expressed as a Laplace transform. Our final 
expression, Eq. (jU3|) . can be interpreted as a Fourier integral, viz (fT3|). which (loosely 

§ In Maldacena's normalization, the numerical prefactor in Eq. © is not consistent with the square 
of the two-point function, Eq. Q. We choose A so that the prefactor becomes 47r^(27r)'^. This 
normaUzation of Eq. ^ was also employed in Refs. |42ll41j . although the distinction from Maldacena's 
A was not pointed out explicitly. Throughout we work in units where the reduced Planck mass, defined 
by Mp = (SttG)^^, is set to unity. If necessary, a finite Planck mass can be restored in formulas such 
as © by dimensional analysis. 
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speaking) can be related to a Laplace integral via a Wick transformation. Despite 
these similarities, the correspondence between the two is complicated because Refs. 
[OHl innj work in a multiple-field picture and calculate a probability density only for the 
isocurvature field '5s', which acquires its non-Gaussianity via a mixing of isocurvature 
and adiabatic modes long after horizon exit. This contrasts with the situation in the 
present paper, where we restrict ourselves to a single-field scenario and compute the 
probability density for the adiabatic mode 7^, which would be orthogonal to 5s in field 
space and whose non-Gaussianity is generated exactly at horizon exit. 

In the older literature it is more common to deal with the density fiuctuation 5p/ p 
measured on comoving slices, as an alternative to the curvature perturbation TZ. For 
slowly varying fields, on large scales, TZ and 5p can be related via the rule (Eq. (25) of 
Ref. [Zq!), 

-H\'5p f3 , 1 Y'^^ ^g^ 



k J p \2 1 +CJ 

which is valid outside the horizon, and to first order in cosmological perturbation theory, 
for a fiuid with equation of state p = up. (One may use the SN formalism to go beyond 
leading order, but to obtain results valid on sub-horizon scales one must use the full 
Einstein equations directly; see, eg., [ZU [Z21-) For fiuctuations on the Hubble scale, 
where k ~ aH, this means |7^| ~ Sp/p, so TZ provides a useful measure of the density 
fiuctuation on such scales. In virtue of this relationship with the density fiuctuation, 
the probability distribution P(e) is an important theoretical tool, especially in studies 
of structure formation. For example, it is the principal object in the Press-Schechter 
formalism jHZI. As a result, there are important reasons why knowledge of the detailed 
form of the probability distribution of e is important, and not merely the approximate 
answer provided by the central limit theorem. 

Firstly, large collapsed objects, such as primordial black holes (PBHs) naturally 
form on the high-e tail of the distribution [7^ I74j . Such large fiuctuations are extremely 
rare. This means that a small change in the probability density for |e| ^ can make a 
large difference in the mass fraction of the universe which collapses into PBHs |75| l76j. 
Thus one may hope to probe it using well-known and extremely stringent constraints on 
PBH formation in the early universe [771 HH] • The corrections calculated in this paper 
are therefore not merely of theoretical interest, but relate directly to observations, and 
have the potential to sharply discriminate between models of infiation. 

Secondly, as described above, although the non-Gaussianities produced by single- 
field, slow-roll infiation are small, this is not mandatory. In models where non- 
Gaussianities are large, it will be very important to account for the effect of non-Gaussian 
fiuctuations on structure formation [021 IHZl HOI IHO] ■ The formalism presented in this 
paper provides a systematic way to obtain such predictions, extending the analysis 
given in Ref. [HZ]. 

The outline of this paper is as follows. In Section |21 we obtain the probability 
measure on the ensemble of possible fiuctuations. This step depends on the correlation 
functions of 7?.. In Section O we discuss the decomposition of 7^ into harmonics. This 
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is a technical step, which is necessary in order to write down a path integral for P(e). 
First, TZ is decomposed in Section 13.11 The path integral measure is written down 
in Section \^.2\ and in Section 13.31 we give a precise specification of e, which measures 
the size of fiuctuations. We distinguish two interesting cases, a "total fiuctuation" e, 
which corresponds to TZ (or approximately Sp/p) smoothed over regions the size of the 
Hubble volume, and the "spectrum" g{k), which describes the contributions to e from 
regions of the primordial power spectrum around the scale described by wavenumber 
k. In Section m we evaluate P(e). We give the calculation for the Gaussian case first, 
in Section 14.11 in order to clearly explain our method with a minimum of technical 
distractions. This is followed in Section ^21 by the same calculation but including non- 
Gaussian corrections which follow from a non-zero three-point function. In Section El 
we calculate P[^(/i;)]. Finally, we state our conclusions in Section IHl 

2. The probability measure on the ensemble of TZ 

Our method is to compute the probability measure P( [R] on the ensemble of realizations 
of the curvature perturbation R(x), which we define to be the value of TZ{t,x) at some 
fixed time t. This probability measure is a natural object in the Schrodinger approach to 
quantum field theory, where the elementary quantity is the wavefunctional ^I/t[f^]) which 
is related to Pt[R] by the usual rule of quantum mechanics, that Pt[R] oc |^i[R]p. Once 
the measure Pt[R] is known, we can directly calculate (for example) Pf(e) by integrating 
over all R that produce fiuctuations of amplitude e. Although the concept of a probability 
measure on R may seem like a rather formal object, the Schrodinger representation || of 
quantum field theory is entirely equivalent to the more familiar formulation in terms 
of a Fock space. Indeed, a similar procedure has been discussed by Ivanov [70], who 
calculated the probability measure on a stochastic metric variable ais(x) which can be 
related to our 7?.(x). Although the approaches are conceptually similar, our method 
is substantially different in detail. In particular, the present calculation is exact in the 
sense that we make no reference to the stochastic approach to infiation, and therefore are 
not obliged to introduce a coarse-graining approximation. Moreover, Ivanov's analysis 
appeared before the complete non-Gaussianity arising from 7^-field interactions around 
the time of horizon crossing had been calculated jlU] , and therefore did not include this 
effect. 

2.1. The generating functional of correlation functions 

The expectation values (7^(xi) ■ ■ ■TZ{^2)) in the vacuum |in) at some fixed time t can 
be expressed in terms of a Schwinger-Keldysh path integral,^ 

(7^(t,Xl)■■■7^(t,x„)) = 

II The Schrodinger representation is briefly discussed in, for example, Refs |81ll82j . A brief introduction 
to inflnitc-dinicnsional probability measures is given in Ref. 1881 . 

^ Henceforth, we use the notation (• • •) to mean expectation values in the in-vacuum, and no longer 
write I in) explicitly where this is unambiguous. 
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y"[d7^_d7^+]^4(*''')=^-(*''')7^(t,Xl) • ■ •7^(t,x„) exp {iS[n+] - iS[n^]) . (9) 

(For details of the Schwinger-Keldysh or "closed time path" formalism, see Refs. 
|T7| UHl EBl 1211 1201 I22I-) in cosmology we are generally interested in TZ evaluated at 
different spatial positions on the same t-slice, so we have set all the t equal in ^. The 
path integral is taken over all fields TZ which begin in a configuration corresponding to 
the vacuum |in) at past infinity, and 5* [7^] is the action for the fluctuation 7^, which is 
given to third order in TZ in Refs. |^ 021 1^ • 

An expression equivalent to Eq. (jHl) can be given in terms of the "equal time" 
generating functional Zt[ri\ at time t, defined by 

z,M = j[dR] y"[d7^_d7^+];f„f*''^)="(")exp Us[n+] - iS[n^] + ijd^x r{-kU^)\ , (lo) 

where rj is some arbitrary source field and S^ is a spatial slice at coordinate time t. This 
generating functional is not the usual one, which would generate correlation function at 
any given set of times, and not the common time t which appears in (J9|l- (|T0|l . The equal- 
time correlation functions (7^(t,xi) ■ ■ ■7?.(t,x„)) are recovered from Zi[ri\ by functional 
differentiation. 



1 A A 

(7^(t,Xl)■■■7^(t,x„))= --—-■■■-— -In ZiM 

i"d?7(xi) dr]{^n) 



(11) 

r;=0 

Up to normalization, this is merely the rule for functional Taylor coefficients, so it is 
straightforward to invert Eq. (jllI) for Zt[ri\. We obtain"*" 



^^ -n p r 

Zt[7]] = expJ2- ■■■ d=^Xl■■■d3a;„r/(xl)■■■r7(x„)(7^(t,Xl)■■■7^(t,x„)), (12) 

Eq. (fTUI) for the generating functional can be rewritten in a suggestive way. We 
define the wavefunctional at time t, "^tiR], as 

n^] = y"[d7^]^i;•'')="w exp {ism . (13) 

This definition is simply the functional generalization of the familiar quantum- 
mechanical wavefunction. It expresses the amplitude for the field TZ{t, x) to have the 
spatial configuration R(x) at time t, given the boundary condition that 7?. started in 
the vacuum state in the far past. In terms of ^l't[R], the generating functional can be 
rewritten 

Zt[r]] = AdR] ^,[R]t^,[R]exp A /"d^x R{^M^)] = \%\R\? « P[R], (14) 

+ This construction is somewhat similar to the dS/CFT calculations outlined in Refs. |4U[ I84L 1551 I86j . 
In these calculations one constructs |^p from an expression of the same form as (|12|l . but expressed in 
terms of an operator O which is the holographic dual of the bulk field TZ. The dS/CFT prescription 
relates the correlators of O reciprocally to those of TZ. In this paper we do not make any use of 
holographic arguments, but the same reciprocal relationship emerges through the identification of Zt[7]\ 
with the Fourier transform of |^p. The normal rules of Fourier transforms show that the transform of 

a Gaussian is another Gaussian with a reciprocal coefficient, e~"^^/^ — > e^'' /^°. 
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where a tilde denotes a (functional) Fourier transform, and f denotes Hermitian 
conjugation. Eq. p4|) implies that Zt[ri\ is the complementary function [HHj for the 
probability distribution Pj[R]) which can formally be obtained by inversion of Zt[ri\. 
Hence, up to an overall normalization. 



Pi[R] oc /[dr/] exp ( -i / d^x R(x)r/(x) J Zt 



W- (15) 



The normalization is not determined by this procedure, but it is irrelevant. We will fix 
the R-independent prefactor, which correctly normalizes the probability distribution, by 
requiring J deP(e) = 1 at the end of the calculation. For this reason, we systematically 
drop all field-independent prefactors in the calculation which follows. 

2.2. The probability density on the ensemble 

So far, all our considerations have been exact, and apply for any quantum field TZit, x). 
For any such field, Eq. |T5|) gives the probability density for a spatial configuration R at 
time t, and implies that to obtain Pt[R] we should need to know all such functions for 
all n-point correlations, and at all spatial positions x. In practice, some simplifications 
occur when TZ is identified as the inflationary curvature perturbation. 

The most important simplification is the possibility of a perturbative evaluation. 
The dominant mode of the CMB fluctuation is constrained to be Gaussian to high 
accuracy, so the corrections to the leading order Gaussian result cannot be large. 
Moreover, since we assume 7?. is the curvature perturbation which is communicated 
to CMB fluctuations, the amplitude of its spectrum is constrained by observation. 
Specifically, in the region of wavenumbers probed by the CORE DMR instrument 
IHH (and other later CMB experiments such as BOOMERANG Hj and WMAP 0), 
the spectrum has amplitude "P^/^ ~ 10~^, whereas the requirement that inflation not 
overproduce PBHs generally requires V^'"^ < 10^^ over the relevant wavenumbers |77j . 
Each higher-order correlation function is suppressed by an increasing number of copies 
of the spectrum, V{k). Provided the amplitude of V is small, it is reasonable to believe 
that we are justified in truncating the exponential in (fTH|l at a given level in n and 
working with a perturbation series in V. The relevant criterion is the smallness of the 
spectrum, rather than the validity of slow-roll. 

However, this simple approach is too naive, because the integrals over t] eventually 
make any given term in the series large, and invalidate simple perturbative arguments 
based on power counting in V. The perturbation series can only be justified a posteriori, 
a point to which we will return in Section 14.21 

We work to first-order in the three-point correlation, 

PM^JM'^tlvMm^l (16) 

where T[ri] and uj[ri; R] are defined by 

(i r cPk d^k cl^k \ 

1 - g y \^^^, ' r/(kl)r/(k2)r^(k3)(7^(t, k,)n{t, k,)n{t, k3))J , (17) 
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and 

.,h; R| = exp (- / ^H^ !?<!^(K(, kOK(t. fe)) - i / jgs .Kk)R(k)) ,(18) 

In this expression, uj will give rise to the Gaussian part of the probability distribution, 
and T is of the form 1 plus a correction. This correction is small when the perturbative 
analysis is valid. Higher-order perturbative corrections in V can be accommodated if 
desired by retaining higher-order terms in the power-series expansion of the exponential 
in ()15|). Therefore our method is not restricted to corrections arising from non- 
Gaussianities described by three-point correlations only, but can account for non- 
Gaussianities which enter at any order in the correlations of TZ, limited only by the 
computational complexity. However, in this paper, we work only with the three-point 
non-Gaussianity. 

We now complete the square in uJt[ri\ R] and make the finite field redefinition 

where the prime ' attached to (7?.(t, k)7?-(t, — k))' indicates that the momentum- 
conservation ^-function is omitted. The measure [drj] is formally invariant under this 
shift, giving j[&q\ = /[dr/], whereas ujt[i]] R] can be split into an R-dependent piece, 
which we call ri[R], and a piece that depends only on f) but not R, 

"2/ ^^^^(kl)^(k2)(7^(t,k07^(t,k2))j ,(20) 

where Tt[R] is a Gaussian in R, 

r,|R] = exp (-1 / d% d-'fe m, kOK«. fa)) j3_^Xk.)^'('.-k.))- ) ■ pi' 

Eq. (fTB|) for the probability density becomes 

P,[R] oc r,[R] J[dv] T, exp (^-^ J ^^0^ ^(ki)r)(k2) (7^(t, k07^(t, k^))) , (22) 

One can easily verify that this is the correct expression, since if we ignore the three-point 
contribution (thus setting T^ = 1), one recovers (after applying a correct normalization) 

AdR] R(ki)R(k2)r4R] = (7^(t,kl)7^(t,k2)). (23) 

The remaining issue is to carry out the fj integrations in Tj. The only terms which 
contribute are those containing an even power of fj, since any odd function integrated 
against e~^ vanishes identically. In the expansion of Yl- //(kj) in terms of 17, there are 
two such terms: one which is quadratic in 17, and one which is independent of fj. These 
are accompanied by linear and cubic terms which do not contribute to Pi[R]- For any 
symmetric kernel K and vectors p, q G R"*, one has the general results 

J[df] exp(^-iy"d"^xd"^y/(x)/(y)K(x,y)') =(detK)-^/% (24) 

y"[d/]/(p)/(q)exp(^-iy"d"^xd"^2//(x)/(y)K(x,y)^=K-l(p,q)(detK)-^/^ (25) 
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These rules allow us to evaluate the r) integrals in Eq. (J22I), giving 

P4R] oc r,[R] (l + Tf\R] + Tf)[R]) , (26) 

where T'^^-' and T*^^-* are defined by 

+ combinations, (28) 

where in T\ we include the possible combinations of the labels {1, 2, 3} which give rise 
to distinct integrands. 

In fact, T\ is negligible. This happens because the three-point function contains a 
momentum-conservation ^-function, 5(ki + k2 + ks), which requires that the vectors k, 
sum to a triangle in momentum space. [For this reason, it is often known as the "triangle 
condition", and we will usually abbreviate it schematically as 5(A).] In combination 
with the (5-function 5(k2 + ks), the effect is to constrain two of the momenta (in this 
example k2 and ka) to be equal and opposite, and the other momentum (in this example, 
ki) to be zero. This is the extreme squeezed limit [101 |HH1 Ell , in which the bispectrum 
reduces to the power spectrum evaluated on a perturbed background, which is sourced 
by the zero-momentum mode. Written explicitly, TJ behaves like 

Tf\R] ~ - / , \, ' ctR(ki)5(ki) + combinations, (29) 

6 J [2-kY 

where we have written linifci^o-^ = ct^i^ ^^^ a is a known finite quantity. In particular, 
Eq. ()29p vanishes, provided R(k) approaches zero as A; ^ 0. This condition is typically 
satisfied, since by construction R(k) should not contain a zero mode. Indeed, any zero 
mode, if present, would constitute part of the zero-momentum background, and not 
a part of the perturbation R. Accordingly, Eqs. (J26|l - (j27jl with T^ = give Pt[R] 
explicitly in terms of the two- and three-point correlation functions. 

2.3. The smoothed curvature perturbation 

The probability density Pt[R] oc (1 + TJ ^[R])rt[R] that we have derived relates to the 
microphysical field 7^(t, x) which appeared in the quantum field theory Lagrangian. A 
given k-mode of this field begins in the vacuum state dX t ^ — oo. At early times, the 
mode is far inside the horizon {k ^ aH). In this ( "subhorizon" ) regime, the k-mode 
cannot explore the curvature of spacetime and is immune to the fact that it is living 
in a de Sitter universe. It behaves like a Minkowski space oscillator. At late times, 
the mode is far outside the horizon {k -C aH). In this ( "superhorizon" ) regime, the 
k-mode asymptotes to a constant amplitude, provided that only one field is dynamically 
relevant during inflation pUl I13j.* If we restrict attention to tree-level diagrams, then 

* Where multiple fields are present, there will typically be an isocurvature perturbation between them: 
hypersurfaces of constant pressure and density will not coincide. Under these circumstances |13| . TZ 
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under reasonable conditions the integrals which define the expectation values of TZ are 
typically pHl |23 dominated by the intermediate ("horizon crossing") regime, where 
7^(k) is exiting the horizon [k ~ aH). As a result, the correlation functions generally 
depend only on the Hubble and slow-roll parameters around the time of horizon exit. 

The simple superhorizon behaviour of TZ means that we can treat the power 
spectrum as constant outside the horizon. As has been described, its value depends 
only on the Hubble parameter and the slow-roll parameters around the time that the 
mode corresponding to k exited the horizon. For this reason, the time t at which we 
evaluate the wavefunctional ^E^J^^]' ^^^ generating functional Zt[ri\ and the probability 
distribution Pt[R] is irrelevant, provided it is taken to be late enough that the curvature 
perturbation on interesting cosmological scales has already been generated and settled 
down to its final value. Indeed, we have implicitly been assuming that t is comoving 
time, so that observers on slices of constant t see no net momentum fiux. Because TZ is 
a gauge invariant, and is constant outside the horizon, our formalism is independent of 
how we choose to label the spatial slices. The evolution of IZ outside the horizon is the 
principal obstacle which would be involved in extending our analysis to the multiple-field 
scenario. 

When calculating the statistics of density fiuctuations on some given lengthscale 
2Ti/kHi OI16 should smooth the perturbation field over wavenumbers larger than ku. To 
take account of this, we introduce a smoothed field R which is related to R via the rule 
R(k) = >V(A;, /cj:/)R(k), where W is some window function. The probabilities we wish to 
calculate and compare to the real universe relate to R rather than R. The exact choice 
of filter W is mostly arbitrary. For the purpose of analytical calculations, it is simplest 
to pick a sharp cutoff in k-space, which removes all modes with k > kn- However, 
this choice has the disadvantage that it is non-local and oscillatory in real space, which 
makes physical interpretations difficult. The most common alternative choices, which 
do not suffer from such drawbacks, are: (i) a Gaussian, or: (ii) the so-called "top hat," 
which is a sharp cutoff in real space. We allow for a completely general choice of C° 
function W, subject to the restriction that W 7^ except at k = 00 and possibly at an 
isolated set of points elsewhere. This restriction is made so that there is a one-to-one 
relationship between R and R. If this were not the case, it would be necessary to coarse- 
grain over classes of microphysical fields R which would give rise to the same smoothed 
field R. 

In addition to this smoothing procedure, the path integral must be regulated 
before carrying out the calculation in the next section. This is achieved by artificially 
compactifying momentum space, so that the range of available wavenumbers is restricted 
to fc < A, where A is an auxiliary hard cutoff. jj At the end of the calculation one 

will evolve. We do not consider the evolving case in this paper, but rather restrict our attention to the 
single-field case where the superhorizon behaviour of TZ is simple. 

jj Note that this procedure does not have anything to do with the regularization of ultraviolet 
divergences. Such divergences are connected to the appearance of loop graphs, which we ignore, and 
in any case are subdominant [23 ^M ■ 
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takes A — i> oo. Some care is necessary in carrying out this compactification. We 
set R = for k > A. In order to maintain continuity at A; = A, we introduce a 1- 
parameter family of functions Wa. These functions are supposed to satisfy the matching 
condition hmA^oo VVa(A;) = W{k), and are subject to the restriction Wa(A) = 0. (These 
conditions could perhaps be relaxed.) The relationship between R and R becomes 

R(k) = 9{A-k)WA{k;kH)R(k). (30) 

where 6{x) = j_ 6{z) dz is the Heaviside function. To minimise unnecessary clutter in 
equations, we frequently suppress the A and kn dependence in W, writing only W{k) 
with the smoothing scale kn and hard cutoff A left implicit. Both of the standard window 
functions approach zero as fc — > cxd, and are compatible with ()3n|l in the A ^ oo limit. 
In this limit, the final result is independent of the exact choice of family VVa(A;, kn)- 

The probability in which we are interested is that of observing a given filtered field 
R. One can express this via the rule (see also |?)7|l91j) 

Pt[R] = AdR] Pt[R]6[R = 9{A - k)WR]. (31) 

3. Harmonic decomposition of the curvature perturbation 

In the previous section, we obtained the probability density for a given smoothed 
spatial configuration of the curvature perturbation. Given this probability density, 
the probability P(C) that the configuration exhibits some characteristic (, such as 
fluctuations of size e or a 'fluctuation spectrum' (to be defined later) of the form g{k), is 
formally obtained integrating over all configurations of R which exhibit the criteria which 
define ( (cf. Ref. |HZ|)- We give a precise specification of these criteria in Section EISl 
below. Before doing so, however, we exploit the compactification of momentum space 
introduced in (jHUj) to define a complete set of partial waves. The smoothed field R can be 
written as a superposition of these partial waves with arbitrary coefficients. Moreover, 
the path integral measure can formally be written as a product of conventional integrals 
over these coefficients [22] • 

In this section we assemble the necessary formulae for the partial-wave 
decomposition. In particular, we shall require Eq. ()32|1 for the decomposition of R, 
Eq. (jH^ for the path integral measure, and Eq. (jUj), which gives a precise specification 
of the characteristics e and [){k). 

3.1. Harmonic expansion of R 

We expand R(k) in harmonics on the unit sphere and along the radial fc = |k| direction, 

oo <? oo 

R(k) = E E E ^"inY^rnie, <P)Mk) ■ (32) 

e=0 m=-i n=l 
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The Y£m,(^!0) are the standard spherical harmonics on the unit 2-sphere (see, eg., Ref. 
|93j). whereas the ipn{k) are any complete, orthogonal set of functions on the interval 
[0,A]. These harmonics should satisfy the following conditionsiff 

(i) ipn{k) -^ smoothly a.s k —>■ 0, so that power is cut off on very large scales, and 
the universe remains asymptotically FRW with the zero-mode a{t) which was used 
when computing the expectation values {TZ- ■ -IZ) left intact; 

(ii) ipn{k) -^ smoothly as fc — » A, so that the resulting R is compatible with Eq. (jHUjl : 

(iii) ipn{k) should have dimension [M~^], in order that Eq. ()32|) is dimensionally correct; 

and 
(iv) the ipn{k) should be orthogonal in the measure J„ dk k^V^^{k)W^'^{k). 

In addition, there is a constraint on the coefficients RT, because R(k) should be real 
in configuration space and therefore must obey the Fourier reality condition R(k)* = 
R(— k), where a star '*' denotes complex conjugation. The R^^ are generically complex, 
so it is useful to separate the real and imaginary parts by writing R^^ = a^i + i^Zi- The 
condition that R is real in configuration space implies 

aJ^Z = i-lY'-'^a^n (33) 

6-^ = i-lY^-'-^'b^n- (34) 

These conditions halve the number of independent coefficients, since the a and b 
coefficients with strictly negative m are related to those with strictly positive m, whereas 
for the m = modes, the b coefficients vanish if i is even, and the a coefficients vanish 
if i is odd. 

Condition (ji|) is made so that the power on large scales is smoothly cut off. In 
the absence of this constraint, R could develop unbounded fluctuations on extremely 
large scales which would renormalize a{t). Therefore, Condition (ji|) can be interpreted 
as a consistency requirement, since the inflationary two- and three-point functions are 
calculated using perturbation theory on an FRW universe with some given a{t), which 
must be recovered asymptotically as |x| -^ oo. It will later be necessary to sharpen this 
condition to include constraints on the behaviour of V{k) near k = 0, beyond the weak 
requirement that a^ = J V{k) d In /c is finite. Condition (jivjl is a technical condition made 
for future convenience. Any other choice of normalization would work just as well, but 
this choice is natural given the /c-dependence in the Gaussian kernel r[R]. Indeed, with 
this condition, the Gaussian prefactor in P(e) will reduce to the exponential of a sum of 

ft When expanding functions on R in terms of polar coordinates, a more familiar expansion is in terms 
of the spherical waves Z£„i\k oc ji{kr)Yijn{S, (f>), where ji is a spherical Bessel function. These waves are 
eigenfunctions of the Laplacian in polar coordinates, viz, V^Z^„i|fe — — fc^Zf„i|fc. An arbitrary function 
on R can be written in terms of the spherical waves, which is equivalent to a Fourier expansion. We do 
not choose the spherical waves as an appropriate complete, orthogonal set of basis functions because we 
do not wish to expand arbitrary functions, but rather functions obeying particular boundary conditions, 
specifically, at fc = 0. The spherical waves for low £ behave improperly at small k for this purpose. 
Moreover, it is not possible to easily impose the boundary condition R(fc) ^- as fc ^ A. 
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squares of the a^^ and 6^^. We stress that in virtue of Condition (jm)), the inner product 

of two ipn{k) in the measure J^ dkk^V~^{k)}V~'^{k) is dimensionless. Condition (jnl) 
has less fundamental significance. It follows from the condition VVa(A) = and the 
artificial compactification of momentum space. However, as in the usual Sturm-Liouville 
theory |M], the precise choice of boundary condition is immaterial when the regulator 
is removed and A — > cxd. Condition (jHl) does not affect the final answer. 

To demonstrate the existence of a suitable set oiipnik), we can adopt the definition 



^ (k) = ^ vjkmk) f k_ 



(35) 



where Ju{z) is the Bessel function of order z/, which is regular at the origin, and a" is its 
nth zero. The order v is arbitrary, except that in order to obey condition (ji]) above, we 
must have k^~'^V{k) -^ as A; — > 0, assuming that W(fc) — > 1 as fc — ;► (as is proper for 
a volume- normalized window function). The ipn{k) obey the orthonormality condition 



dk 



-i)n{k)i)m{k) 



and 5r, 



,0 V{k)W^{k) 

is the Kronecker delta. The completeness relation can be written 



Kk - kQ)\ke[o,K] 



l5 



V{k^)W\k, 



^ipnik)i;niko), 



(36) 



(37) 



where the range of k is restricted to the compact interval [0, A]. 

Although we have given an explicit form for the ipn in order to demonstrate 
existence, the argument does not depend in detail on Eq. (jH^ . The only important 
properties are Eqs. (|H ^ -(|H7 |) . which follow from Condition (jlvj). 



3.2. The path integral measure 



k^O 



Since any real, C° function R obeying the boundary conditions R(k) ^^-^ and 
R(k) ^—>- can be expanded in the form (j32j) . one can formally integrate over all 
such R by integrating over the coefficients R^^. This prescription has been widely used 
for obtaining explicit results from path integral calculations. (For a textbook treatment, 
see Ref. 1^51 •) In the present case, it must be remembered that one should include in 
the integral only those R(x) which are real, and can therefore correspond to a physical 
curvature perturbation in the universe. Since the Y^^ sue complex, this means that 
instead of integrating unrestrictedly over the R^^, the reahty conditions (I^Hj) must be 
respected. A simple way to achieve this is to integrate only over those a^^ or &T with 
m > 0. The m = modes must be treated specially since the a and b coefficients vanish 
for odd and even £, respectively. 

The integral over real R can now be written 



[dR] 



R 



nnn" 

£=0 m=l n=l 



da 



■e\n 






r=0 s=l 



da 



d/)0 



{31 
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where the subscript R on the integral indicates schematically that only real R(x) are 
included. The constants /i and fx account for the Jacobian determinant which arises 
in writing /[dR] in terms of the harmonic coefficients RT. Their precise form is of no 
importance in the present calculation. 

As we noticed above, the detailed form of the measure (jHH|) is not absolutely 
necessary for our argument. The important aspect is that each a or 6 integral (where 
m > 0) can be carried out independently. For this purpose it is sufficient that the 
spectrum of partial waves is discrete, which follows from the (artificial) compactness of 
momentum space. However, although it is necessary to adopt some regulator in order 
to write the path integral measure in a concrete form such as (jHH|) . we expect the answer 
to be independent of the specific regulator which is chosen. In the present context, this 
means that our final expressions should not depend on A, so that the passage to the 
A — > oo limit becomes trivial. 

3.3. The total fluctuation e and the spectrum g{k) 

There are at least two useful ways in which one might attempt to measure the strength 
of fiuctuations in R. The first is the total smoothed fluctuation at a given point x = xq. 
By a suitable choice of coordinates, we can always arrange that xq is the origin, so the 
condition is R(0) = e. When R is smoothed on scales of order the horizon size this 
gives a measure of the fluctuation in each Hubble volume, since distances of less than 
a horizon size no longer have any meaning. For example, Shibata & Sasaki jHSj have 
proposed that e defined in this way represents a useful criterion for the formation of 
PBHs, with formation occuring whenever e exceeds a threshold value eth of order unity 
l^j . This measure of the fiuctuation is non-local in momentum space. Making use of the 
relation Jdf2(6',0) Yimid,4>) = V^^efi^mfl ^oi the homogeneous mode of the spherical 
harmonics, one can characterize e as 

j^.me'-^^o = ^,JdkeY. <nUk) = e. (39) 

On the other hand, one might be interested in contributions to the total smoothed 
fiuctuation in each Hubble volume which arise from features in the spectrum near some 
characteristic scale of wavenumber k. For this reason, we consider a second possible 
measure of fiuctuations, which we call the fluctuation spectrum, and which is defined 
by the condition g{k) = dR(0)/dlnA;. (Thus, the total smoothed fiuctuation can be 
obtained by integrating its spectrum according to the usual rule, viz, e = J g{k) dlnfc.) 
This condition is local in k-space. Differentiating (J39|l . one can characterize g{k) as 



^(k) = -^J2^o\nk'Mk)- (40) 

This is a functional constraint. 

We will calculate the statistics of both the total fiuctuation e and its spectrum g{k). 
In each case, the calculation is easily adapted to other observables which are non-local 
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or local in momentum space, respectively. Indeed, both the non-local e and the local 
Q{k) are members of a large class of observables, which we can collectively denote C,, 
and which all share nearly-Gaussian statistics. Specifically, Eqs. ()H9|1 and ^^ can be 
written in a unified manner in the form 

J]a°|„E„(fc) = ^C(fc), (41) 

— : v47r 

n=l * 

where the S„ are defined by 

Jp dk k'^ipri(k) total fluctuation, C = e; , . 

k^ipn{k) fluctuation spectrum,^ = g{k). 

Note that in the case of e, the E„ are independent of k. Any characteristic which can be 
put in this form, with a coupling only to the real zero-modes Ogi of R, will necessarily 
develop nearly-Gaussian {i.e., weakly non-Gaussian) statistics. More general choices of 
characteristic are possible, which cannot be cast in the form (J4H) . For example, one can 
consider characteristics which depend non-linearly on the a^.. Such characteristics 
will generally lead to strongly non-Gaussian probabilities. The Gaussianity of the 
final probability distribution depends on the geometry of the constraint surface in an 
analogous way to the decoupling of the Fadeev- Popov ghost fields in gauge field theory 
PH] . These non-Gaussian choices of characteristic can also be handled by generalizing 
our technique, but we do not consider them here. 

4. The probability density function for e 

We first calculate the probability density for the non-local constraint e, given by Eq. PUJ). 
The expression is 



P(e) oc / [dR] P[R]6 
Jn 



, V47r 

n=l * 



(43) 



To obtain this density, one treats e as a collective coordinate parametrizing part of 
R. The remaining degrees of freedom, which are orthogonal to e, are denoted R^. 
Therefore the functional measure can be broken into [dR] oc [dR^] de. After integrating 
the functional density P[R] [dR] over R^, the quantity which is left is the probability 
density P(e) de. In this case, the integration over the orthogonal degrees of freedom R-*- 
is accomplished via the 5-function, which filters out only those members of the ensemble 
which satisfy Eq ()39|) . We emphasize that this is a conventional (5-function, not a 6- 
functional. There is no need to take account of a Fadeev-Popov type factor because 
the Jacobian associated with the constraint (PT|) is field- independent, in virtue of the 
linearity of (jHnjl in a^,^. 

4.1. The Gaussian case 

We first give the calculation in the approximation that only the two-point function 
is retained. In this approximation, the probability distribution of e will turn out to 
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be purely Gaussian, which allows us to develop our method without the distractions 
introduced by including non-Gaussian effects. 

If all correlation functions of order three and higher are set to zero, then P[R] oc 
r[R]. Using (0), one can write 

1 /■ ,^ /".o ,, k^ 1 



r[R] =exp(--fdnfk^dk 



{2n)^27i^V{k)W^{k) 



X E E K\nXt'^^^^^^^^^Krn^O,<l^)^nAk)^nM]- (44) 

£i ,mi ,n-i £2 ,^2 ,112 

The harmonics Y^^ and ipn integrate out of this expression entirely, using the 
orthonormality relation ()36|) and the spherical harmonic completeness relation 
J dilY£-^mjY^2m2 ~ ^iii2^rnim2- Morcovcr, after rewriting the a and b coefficients with 
m < in terms of the m > coefficients, we obtain 

(^oofoo ^OOOO \ 

'^ '' e=0 m=l n=l ^ ' e=0 n=l ^ 

£ even 

The 5-function in P3|) constrains one of the a^, in terms of e and the other 
coefficients. It is possible to evaluate P(e) by integrating out the 5- function immediately. 
This would involve solving the constraint for Oqiq (for example) and replacing it in the 
integrand with its expression in terms of the other Qq. and e. However, this does not 
turn out to be a convenient procedure, for the same reasons that one encounters when 
gauge-fixing in field theory. Instead, we introduce the Fourier representation of the 
(5-function, 



Pfeloc 



f [dR] f d^r[R]exp liz 
Jn i-oo V 



E«o|.S, 






(271' 

' e 



V 47r 

.n=l 



(46) 



where the functional measure is understood to be Eq. (jHSj) . The final answer is obtained 
by integrating out z together with all of the a and b coefficients. In order to achieve 
this, it is necessary to decouple Oqi^, z and e from each other by successively completing 
the square in Oqiq and z. Working with Oqiq first, we find 



(^ ^ oo oo > 

~4^TW E I^Ol^r + i^ E «0|nSn 
^ ' n=l 71=1 y 

1 1 

'47r2(27r) 



^^P I -Z^T^ E(«o|n - 2n\2nfizi:^r - (2vr)V^.2s2 , (47) 

n=l / 

where we have introduced a function S^, defined by S^ = J2'^=i ^n- ■'■^ ^^e final 
probability distribution, S^ will turn out to be the variance of e. From Eq. (j47p . it 
is clear that making the transformation a^, ^-^ a^, + 2TT'^{27i)^izT,n suffices to decouple 
Oqi^ from z. The measure, Eq. (jSHI), is formally invariant under this transformation. 
Exactly the same procedure can now be applied to z and e, giving 



exp (-{2n)Vz^J:' - ^^iez) = exp -(27r)37r2s2 (z + '-^—] 

V v^ ; ^ ^ ^ ^ 27r2v^S2; 



> 2 2 



2S2 



.(4^ 
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As before, the finite shift z ^^ z — ie/27r^V47rS^ leaves the measure intact and decouples 
z and e. The a, h and z integrals can be done independently, but since they do not 
involve e they contribute only an irrelevant normalization to P(e). Thus, we obtain 
Gaussian statistics for e, 



^2 



P(e) (X exp ( -^ ) . (49) 



It remains to evaluate the variance S . In the present case, we have S^ 





Jq dk k^ipn{k)- From the completeness relation Eq. (jHTj) . it follows that 



J] klMko)k'Mk) = ^'^(y(^°) ^(fc - ko). (50) 

n " 

E^ is now obtained by integrating term-by-term under the summation. The result 
coincides with the smoothed conventional variance (compare Eq. ^). 

.A 

T.l{kH)= d\nkW\k]kH)V{k). (51) 

Jo 

Thus, as expected, Eq. fH^ reproduces the Gaussian distribution Q which was derived 

on the basis of the central limit theorem, with the proviso that the parameters (such 

as E^) describing the distribution of e are associated with the smoothed field R rather 

than the microphysical field R. Note that E^ is therefore implicitly a function of scale, 

with the scale dependence entering through the window function. 

In particular, it was only necessary to use the completeness relation to obtain this 

result, which follows from Condition (jrvjl in Section ITTl 

4-2. The non-Gaussian case 

The non-Gaussian case is a reasonably straightforward extension of the calculation 
described in the preceding section, with the term T*^°^ in Eq. (j^Bj) (which was dropped in 
the previous section) included. However, some calculations become algebraically long, 
and there are subtleties connected to the appearance of the bispectrum. 

The inclusion of T^^-* corrects the pure Gaussian statistics by a quantity proportional 
to the three-point function, iJZTZTZ), which is given in Eq. (jHl). This correction is written 
in terms of the representative spectrum P^, which describes when the slow-roll prefactor, 
given by the amplitude of the spectrum, should be evaluated ^01- For modes which cross 
the horizon almost simultaneously, /ci ~ A;2 ~ k^, this prefactor should be V^ = V{kY, 
where k is the common magnitude of the fcj. In the alternative case where one k-mode 
crosses appreciably before the other two, V^ should be roughly given by 

p2 = V {max k.i)V {rain ki). (52) 

Since the difference between this expression and V{k)^ when all k are of the same 
magnitude is very small, it is reasonable to adopt Eq. (jHSl) as our definition of P^. We 
stress that this prescription relies on the conservation of TZ outside the horizon jHl], and 
therefore would become more complicated if extended to a multiple field scenario. 
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With this parametrization, the probabihty measure on the ensemble is obtained by 
combining (IKll. (OH|l . (I?fl) and ipi. 

PfRi v\^^(^ 1 f d^hd^k^d^h ^,^^ PU R(kOR(k,)R(k3) ^ 

This expression should be integrated with the constraint (jH^ and measure ()38|) to obtain 
the probability P(e). At first this appears to lead to an undesirable consequence, since 
the integral of any odd function of R against r[R] must be zero. It may therefore seem 
as if the non-Gaussian corrections we are trying to obtain will evaluate to zero. This 
would certainly be correct if the integral were unconstrained. However, the presence 
of the constraint (5-function means that the shifts of a^. and z which are necessary to 
decouple the integration variables give rise to a non-vanishing correction. 

The finite shift necessary to decouple Oq, and z is not changed by the presence 
of non-Gaussian corrections, since it only depends on the argument of the exponential 
term. This is the same in the Gaussian and non-Gaussian cases. After making this 
shift, which again leaves the measure invariant, the integration becomes 

P(e) oc [ [dR] r dz r[R] exp ( -{2nfn^J:h^ - ^^i^e^ (1 - Jo - J2), (54) 

Jr J-oo V v47r / 



where Jq is given by 



I St.St.,%'^HA)^ E ''^'^".^..^ 



712 ^ns 



W{h)Wik2)Wik3 



ni,n2,n3 

and J2 is 

d^kid^k2d^ks P^A y^ ^p >;p 

6(27r)3v/4^ ^ ^U^nh)^,^ ,^ 

^ ^ * "1 ^2,^2,712 ^3 .ms, "2 

+ [1 ^ 21 + II ^ 31. (55) 

The symbol |1 ^ 2l represents the expression in square brackets [■ ■ ■] with the label 
1 exchanged with the label 2, and similarly for |1 ^ 31 . The range of the m2, m^ 
summations is from —£2 to £2 and —£3 to £3, respectively. In addition, the shift of Oq. 
generates other terms linear and cubic in the RT , but these terms do not contribute to 
P(e) and we have omitted them from ()54p . 

After shifting z to decouple z and e, the integrals Jq and J2 develop terms 
proportional to z^, z, z^ and z^. Of these, only the z^ and ^^ survive the final z 
integration. Correspondingly, we suppress terms linear and cubic in z from the following 
expressions. The integral Jq becomes 

Jo = / d^fci d^k2 d^/cg ^ ' ' TTT^^ - 3^^ ^(^ 



ni,n2,n3 



3(47r)2 Vl67r5S6 S^ ^ U^'Pik^) 
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whereas J2 simplifies to 



J, 



E E E 



487r3(27r)3 S2 

^"1 VV(^ ^^1"^ ^^1"^ '^"^^ ''^'^ '^"^^ ''^'^ W(fc2)W(fc3) 

+ [1 ^ 21 + |1 ^ 31, (57) 

and the ?7i summations are still over the entire range, —£2 < ^^2 < ^2 (and similarly for 
7713). Thus Jo contains corrections proportional to e and e^, whereas J2 only contains 
corrections proportional to e. 

The a, h and 2; integrations can now be performed, after the integrand has been 
written entirely in terms of the a^^ and Vm^ with ?ti > 0. There are no a or 6 integrations 
in Jq. In J2, there are no z integrations, but the a and h integrations involved in the 
product R^? R™? fix £2 = -^3, ''^2 = ""^3 and ^2 = ^3. One then uses the spherical 
harmonic completeness relation, 
00 i 

Y, E YU0l^<Pl)yLi^2,<p2) = <^(01 - 02)'5(COS^1 - COS ^2) (58) 

£=0 m=-l 

and the equivalent relationship for the ■j/'-harnionics, Eq. ()37|) . to obtain 
d3fcid3fc2d3fc3 e ^, . , V^A V{k2)W\k2) 



J- 



+ [1 ^ 21 + II ^ 31. (59) 



The terms with 1 exchanged with 2 and 3 do not generate quantitatively different 
integrands and can be absorbed into an overall factor of 3. 
Jo involves only z integrations. It can be written 



967r2 V^e ^V ^ Ili^l^i)"^!^*) 

X ^ S„jE„2S„3?/'„^(/i;i)?/'„2(A;2)V^„3(A;3). (60) 

ni,n2,n3 

To simplify these expressions further, it is necessary to obtain the value of the sum 
E^=iS„^n(A;). 

Reasoning as before from the completeness relation Eq. ()37|1 . it follows that 

f;s„,„,*)^^«^. (61) 

n=l 

From this, it is straightforward to show that Jo behaves like 

d^/ci d^/c2d^A;3 ^z ^ N'fS2 ^( ^ o <= 



where T? is the smoothed variance, Eq. (j^Tj) . On the other hand J2 becomes 

/• d^fcl d3fc2 d^fca 6 ^, ^ , P2 ^^,^,^^^ ^(k2 + k3) 
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After integrating out ka and the angular part of ki and k2, this is the same as 

J. = 2.JAk, kljik, i(tO^>V(*=0^ l.m/^. (64) 

where we have used the fact that ki is constrained to zero by the ^-function to evaluate 
the bispectrum A in the 'squeezed' limit where one of the momenta goes to zero 
|4()| ISni IHHj- In this limit, minfcj = ki and maxfcj = ^2 = ^3, so it is possible to 
expand P^ unambiguously. Moreover, limfc^^o-^ = ct^i i^ proportional to fcfi so J2 = 
if V{k)/k^ — > as A; — s> 0. This is the sharper condition on how strongly large-scale 
power is suppressed which was anticipated in Section ITTl It requires that V{k) cuts off 
on long lengthscales faster than k^ . If this does not occur, then the integral diverges. 
(There is a marginal case where V{k)/k^ tends to a finite limit as k approaches zero. 
We assume that this case is not physically relevant.) 

The J2 integral contains a (5-function 6(k2 + ka). It can therefore be interpreted 
as counting contributions to the bispectrum which come from a correlation between 
the modes k2 and ks, in a background created by ki, which exited the horizon in the 
asymptotic past. As we have already argued, modes of this sort are included in the 
FRW background around which we perturbed to obtain the correlation functions of TZ, 
so we can anticipate that its contribution should be zero, as the above analysis shows 
explicitly. In this interpretation, the condition V{k)/k^ -^ as A; — > is the condition 
that the perturbation does not destroy the FRW background. Indeed, fluctuations on 
very large scales in effect describe transitions from one FRW world to another via a shift 
in the zero-momentum modes of the background metric. In this case, there is only one 
such mode, which is the scale factor a{t). These transitions are rather like changing the 
vacuum state in a quantum field theory. As a result, fluctuations of a large volume of 
the universe between one FRW state and another are strongly suppressed. 

For the case of fluctuations on the scale of a Hubble volume, therefore, the 
probability distribution should be written 

1 
V2^S 

where we have used the fact that the corrections are odd in e, and therefore do not 
contribute to the overall normalization of P(e). The (dimensionless) coefficient J is 

This expression is remarkably simple. Indeed, although the explicit expression (jMj) is 
preferable for calculation, it can be recast directly as the integrated bispectrum with 
respect to W: 

•^ = 48(27r)3(27r2)3 / "^'^^ "^'^^ "^'^^ {n{k,)nik,)ni\.,))W ik,)W ik,)W (k;) (67) 

As a consistency check, we note that the expectation of e, defined by E(e) = 
J eP(e) de satisfies E(e) = 0. This is certainly necessary, since the universe must contain 
as many underdense regions as overdense ones, but is a non-trivial restriction, since 



P(e) 



e^ „ e 



l-(^-3^1^ 



e^ 



exp ( -^ ) , (65) 
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both the e and e^ corrections to P(e) do not separately average to zero. The particular 
combination of coefficients in (jU^ is the unique correction [up to O(e^), containing only 
odd powers of e] which maintains E(e) = 0. 

Finally, we note that Eqs. (j^ -(|U ^ do not explicitly involve the cut-off A, except as 
a limit of integration and in quantities such as T\ and Wa which possess a well-defined, 
finite limit at large A. As a result, there is no obstruction to taking the A — »• cxd limit 
to remove the regulator entirely. 

4-3. When is perturbation theory valid? 

It is known from explicit calculation that the bispectrum is of order V^ multiplied 
by a small quantity, /nl, which is predicted to be small when slow-roll is valid. It 
is therefore reasonable to suppose that whenever the window functions W are peaked 
around some probe wavenumber fc^,, one has the approximate relations (on order of 
magnitude) T? ~ Vi, and J ~ P^, where P* represents the spectrum evaluated at 
k = k^. Since the e^ correction dominates for e > -\/3S, this means that for e not too 
large, e -C P^ , the perturbative correction we have calculated will be small. As e 
increases, so that e ^ V^ , perturbation theory breaks down and the power series in e 
needs resummation. In any case, at such high values of e, the calculation described above 
ought to be supplemented by new physics which can be expected to become important 
at high energy density. The details of these corrections presumably do not matter too 
much, because at any finite order, the fast-decaying exponential piece suppresses any 
contributions from large values of e. 

At some value of e, corrections coming from the trispectrum can be expected to 
become comparable to those coming from the bispectrum that we have computed. Let us 
assume that for intermediate values of e, such that S ^ e ^ V^ , a contribution of the 
form e^S~^ Jtri dominates the correction from the trispectrum, where Jtri is the integrated 
trispectrum with respect to W. The trispectrum contribution will be subdominant 
provided 

^,J>^,Jtru (68) 

which is true whenever e < TP'J/J^^i. The trispectrum from slow- roll inflation has not 
yet been explicitly computed, but it is expected to be proportional to three powers of 
V. Therefore, TP'J^Vl and Jtri ~ 'Pi can be expected to be of roughly equal orders of 
magnitude, up to a numerical coefficient which can estimated to be J/ Jtri — IStt ~ 50. 
(This number comes from 4!/3! = 4 and a factor 47r from the extra angular integrations 
which arise when integrating the trispectrum.) Thus for values of e of order unity, it 
is reasonable to expect the correction from the bispectrum to dominate the correction 
from the trispectrum and higher correlation functions. However, we caution that if the 
trispectrum is anomalously large, or for much larger values of e, higher-order terms in 
the perturbation series will become relevant. 
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5. The probability density function for Q{k) 

The probabihty density function for (){k) can be obtained by a reasonably 
straightforward modification of the above argument, taking account of the fact that 
the constraint, Eq. P0|l . is now a functional constraint. This means that when splitting 
the functional measure [dR] into a product of [d^(A;)] and the orthogonal degrees of 
freedom [dR-*-], the result after integrating out the R^ coordinates gives a functional 
probability density in [df)(A;)]. In particular, the 5-function is now represented as 

j[dz] exp iy"dA; z{k) \f2<\nk'Mk) - ^^(^) 1 • (69) 



,n=l 

In order to carry out this calculation, we write z{k) formally as 

The integration measure ^[dz] becomes HnA^/-^ ^^^^ where, as before, /i is a field- 
independent Jacobian representing the change of variables from z{k) i— »• Zn- Its value is 
not relevant to the present calculation. In addition, we introduce a set of coefficients Qn 
to describe ^(/c), 

'-^ = Y.9nUk). in) 

n=l 

The Qn can be calculated using the rule Qn = /q dA; k'^V'^ {k)'W~'^ (k) Q{k)'ilJn{k) ■ Note 
that in order to do so, we have made the implicit assumption that Q{k)/k^ — ^ as 
fc — i> 0, in order that (ffT| is compatible with the boundary conditions of the ipnik)- We 
again see the suppression of power in modes with low k. 
With these choices, the constraint 5-function becomes 



Y\fi dzn exp 



oo 



1 V v47r 

m=l ^ * 



(72) 



As opposed to the nonlocal case of e, where a single extra integration over z coupled 
to e, we now have a situation where a countably infinite tower of integrations over Zn 
couple to the the coefficients Qn- In all other respects, however, this calculation is now 
much the same as the nonlocal one, and can be carried out in the same way. The shift 
of variables necessary to decouple Oq, and z„ is 

«o|n^«o|n + i27r'(27r)'2;„; (73) 

and the shift necessary to decouple the z„ and Qn is 

Zn^Zn- . 74 

27r2V47r 



When only the two-point function is included, we obtain a Gaussian in the ^„, 

(75) 



P[^(/c)] ocexp {-i^^qI\ 
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The sum over the Qn can be carried out using the completeness and orthogonahty 
relations for the ipn{k), and Eq. (f7T| . 

Using this expression, and integrating over all g{k) which give rise to a fluctuation of 
strength e, one recovers the Gaussian probability proflle Eq. (jl^ with variance given 
by Eq. (jSlj) . This serves as a consistency check of (f7H|l and pHj) . 

When the non-Gaussian correction T^"-' is included, one again generates a 
probability density of the form 

P[Qik)] ^ {I - Ko - K2) exp (-^J2 ^n) , (77) 

where K2 is of the same form as J2, and therefore vanishes for the same reasons; and 
Kq has the form 

The flrst term contains a 5-function which squeezes fci into the asymptotic past. It 
formally vanishes in virtue of our assumption about the behaviour of g{k) near k = 0, 
which is implicit in Eq. (t7T| . As a result, the total probability density for the fluctuation 
spectrum can be written 

P[,(A:)] ex (1 -K)exp (-^/dlnfc ^(|^) , (79) 

where K is given by 

f d^kid^k2d^k3 ^, ^2 -pr Q{h) 

As before, one can show that this expression is consistent with Eqs. (jF)K|l - (pT!)|l by 
integrating over all g{k) which reproduce a total fluctuation of size e, after dropping 
another term which is squeezed into the asymptotic past owing to the presence of a 
5-function. This is a non-trivial consistency check of ()79|) - (jHm) . 

As in the local case, Eqs. (fT^ - ljHUj) are entirely independent of A (except as a limit 
of integration), so the regulator can be freely removed by setting A = 00. 

6. Conclusions 

In this paper we have obtained the connexion between the n-point correlation 
functions of the primordial curvature perturbation, evaluated at some time t, such 
as (7?,(ki) ■ ■ ■7^(k„)), and the probability distribution of fluctuations in the spatial 
configuration of TZ. We have obtained an explicit expression for the probability of 
a fluctuation of "size" e when TZ is smoothed over regions of order the horizon size. 
This is a probability density in the conventional sense. In addition, we have obtained 
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an expression for the probability that e has a spectrum of the form Q{k), that is, 
j dink Q{k) = e. (The mapping g{k) ^^ e is many-to-one.) This is a functional 
probability density, and can potentially be used to identify features in the fluctuation 
spectrum near some specific scale, say of wavenumber k ~ A;^. Our result is independent 
of statistical reasoning based on the central limit theorem and provides a direct route to 
incorporate non-Gaussian information from the vertices of the effective quantum field 
theory of the inflaton into theories of structure formation. 

Both these probabilities are Gaussian in the limit where 7?. only possesses a two- 
point connected correlation function. If there are higher-order connected correlation 
functions, then 7?. exhibits deviations from Gaussian statistics, which we have explicitly 
calculated using recent determinations of the inflationary three-point function during 
an epoch of slow-roll inflation. Our method can be extended to incorporate corrections 
from higher connected n-point functions to any finite order in n. We have not computed 
these higher corrections, since we anticipate that their contribution is subdominant to 
the three-point correction (which is already small), and in any case the relevant 4- and 
higher n-point correlation functions are not yet known explicitly. 

Our argument is based on a formal decomposition of the spatial configuration of 
the curvature perturbation in k-space into spherical harmonics, together with harmonics 
along the radial k direction. However, we have emphasized that our results do not 
depend on the details of this construction, but require only a minimal set of assumptions 
or conditions. These assumptions are: (A) that the power spectrum V{k) goes to zero 
sufficiently fast on large scales, specifically so that V{k)/k^ — ;► as A; — ;► 0. (In addition, 
in the case of the fluctuation spectrum, we also require g{k)/k^ — > as /c — * 0.) Such a 
condition is certainly consistent with our understanding of large scale structure in the 
universe, and within the perturbative approach we are using, we have argued that in fact 
it describes a self-consistency condition which prevents perturbative fluctuations from 
destroying the background FRW spacetime. In addition, we require a second condition 
(B) that the spatial conflguration R can be smoothed R via a window function W to 
obtain a conflguration for which R ^ as A; ^ oo, for which it is fair to compare R to 
the primordial power spectrum. 

In addition to these fundamental assumptions, which relate to the behaviour of 
real physical quantities, a large part of the calculation relied on an auxiliary technical 
construction. This construction is based on an artiflcal compactiflcation of momentum 
space, implemented by a hard cutoff A. There is an associated boundary condition 
on R at A; = A which discretizes the harmonics (partial waves) in k. However, in 
both the non-local (total fluctuation, e) and local (fluctuation spectrum, g{k)) cases, 
the flnal probability density is independent of both the details of the partial wave 
construction and A (except as a limit of integration). It is also independent of the 
choice of the family of window functions >Va(A;; A;^;^), and depends only on the limit 
limA->oo VVa(A;; Ajji^-) = >V(A;;A;//). Therefore the regulator can be removed by sending 
A ^ oo. Moreover, the boundary condition at A; = A becomes irrelevant in this limit, 
which is a familiar result from the theory of Sturm-Liouville operators. As a consistency 
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check, one can integrate P[^(/i;)] with the condition J din A; Q{k) = e in order to obtain 
P(e). 

We conclude by relating our result to earlier work by Ivanov [ZEj and Bullock 
& Primack [7S], who studied the effect of non-Gaussianities from inflation on PBH 
formation, but obtained (naively) opposite results. In Ref. jZHI, it was found that non- 
Gaussianities reduce the probability for PBH formation, whereas in Ref. [7S], using a 
different formalism, the formation probability was found to be enhanced. Consider the 
probability density for the total fluctuation, e. The non-Gaussian correction conserves 
probability, so it can either act by moving probability from the tails into the central 
region, or by moving probability from the central region into the tails. (This effect was 
also observed in Refs. [HHlinSl-) In the former case, the probability of large fluctuations — 
and hence the fraction of the universe going into large collapsed objects — is diminished; 
in the latter case, this situation is reversed. Exactly which occurs depends on the sign 
of the coefficient J . When J > 0, the volume of the central region is increased at the 
expense of the tails. When J < 0, the volume of probability in the tails is increased at 
the expense of the central region. The crossover from central region to tails occurs at a 
threshold amplitude e* = a/3S. A priori, it would appear that J can have either sign, 
depending on the model of inflation under consideration. This goes some way towards 
reconciling the apparently divergent results of Refs. [731 I7H]. 
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